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ABSTRACT 

Recent spectropolarimetric observations of the classical T Tauri star BP Tau and anal- 
ysis of its surface magnetic field have shown that the magnetic field can be approx- 
imated as a superposition of slightly tilted dipole and octupole mom ents with re- 
specti ve strengths of the polar magnetic fields of 1.2 kG and 1.6 kG ( |Donati et al.| 
2008 , hereafter D08). We adopt the measured properties of BP Tau and model the 
disc accretion onto the star by performing global three-dimensional magnetohydrody- 
namic simulations. We observed in simulations that the disc is disrupted by the dipole 
component and matter flows towards the star in two funnel streams which form two 
accretion spots below the dipole magnetic poles. The octupolar component becomes 
dynamically important very close to the star and it redirects the matter flow to higher 
latitudes and changes the distribution and shape of the accretion spots. The spots are 
meridionally elongated and are located at higher latitudes, compared with the pure 
dipole case, where crescent-shaped, latitudinally elongated spots form at lower lati- 
tudes. The position and shape of the spots are in good agreement with observations. 

The disk-magnetosphere interaction leads to the inflation of the field lines and to 
the formation of magnetic towers above and below the disk. The magnetic field of BP 
Tau is close to the potential inside the magnetospheric surface, where magnetic stress 
dominates over the matter stress. However, it strongly deviates from the potential at 
larger distances from the star. 

A series of simulation runs were performed at different accretion rates. In one of 
them, the disk is truncated at r w (6 - 7)R* which is close to the corotation radius, 
R C or ~ 7.5ii*. However, the accretion rate, 1.4 x lO _9 M yr _1 , is lower than that 
obtained from most of the observations. In a sample model with a higher accretion 
rate 8.5 x lO _9 M yr _1 , the disk is truncated at r w 3.6R+, but such a state can not 
be a typical state for the slowly rotating BP Tau if it is in the rotational equilibrium. 
However, torque acting on the star is also small: it is about an order of magnitude 
lower than that which is required for the rotational equilibrium. We suggest that a 
star could lose most of its angular momentum at earlier stages of its evolution. 

Keywords: accretion, accretion discs - magnetic fields - MHD - stars: magnetic fields. 



1 INTRODUCTION 



Accretion-powered Classical T Tauri stars (CTTSs) are young 
low-mass stars which often show signs of a strong magnetic 
field (e.g.,|Basri, Marcy & Valenti|1992[|Johns-Krull, Valenti &| 
Koresko 1999 ) which is expected to have a complex structure 
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(e.g., Johns-Krull 2007). The Zeeman-Doppler imaging tech- 
nique has proven very successful in obtaining surface mag- 
netic maps for many stars, and the external magnetic fields 
of the stars have been reconstructed from these maps under 
the potential appro ximation ([Donati & Col lier Cameron 1997 



Donati et al.|1999| |Jardine et al.||2002||2006[ |Gregory etaL] 
2010). The magnetic field plays a crucial role in disc accretion 
by disrupting the inner regions of the disc and channeling the 
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matter onto the star, and hence it is important to know the 
magnetic field configurations in magnetized stars. 

D08 recently observed the CTTS BP Tau with the ES- 
PaDOnS and NARVAL spectropolarimeters and reconstructed 
the surface magnetic field from the observations. They have 
shown that the magnetic field of BP Tau can be approximated 
by a combination of dipole and octupole components of 1.2 
kG and 1.6 kG, which are slightly (but differently) tilted about 
the rotational axis. 

D07 analyzed the distribution of the accretion spots on 
the stellar surface and found spots at high latitudes, which 
cover about 8 per cent of the stellar surface (D08) . 

Further, D08 extrapolated the surface magnetic field to 
larger distances using the potential approximation, (i.e., as- 
suming that there are no currents outside the star and hence 
the external matter does not influence the initial configura- 
tion of the field) and estimated the distance at which the 
disk should be disrupted by the magnetosphere so that the 
matter flowing towards the star in funnel streams produces 
the high-latitude spots. They concluded that this distance 
should be quite large, r > 4i?*. However, this problem re- 
quires more complete analysis based on the MHD approach, 
where external currents can be taken into account, and the 
matter flow around the magnetosphere can be investigated 
self-consistently taking into account interaction of the exter- 
nal plasma with the magnetic field. 

In this paper, we investigate this problem using global 
three-dimensional MHD simulations. We solve the 3D MHD 
equations numerically in our simulation model to investigate 
the structure of the external magnetic field, accretion flows 
and location of accretion spots. 

In our previous work, we have performed global 3D 
simulations of accretion onto stars with misaligned dipole 
fields ( [Romanova et al.||2003U2004a[ |Kulkarni & Romanova| 
2005 J, and aligned or misaligned dipole plus quadrupole 
fields ( [Long, Romanova & Lovelace||2007| |2008| ). Recently, 
we were able to extend our method and to build a numer- 
ical model for stars with an octupolar component. The gen- 
eral properties of the model have been described in detail in 
Long, Roman ova & Lamb] p010 ) and the model was applied 
to another CTTS, V2129 Oph, with a strong octupole field in 
Romano va et al.| ( [201 0|) which has been compared with ob- 
servations of V2 1 29 dDonati et al.|2007| |Donati et al.|2010| ) . 

In this paper, we apply our 3D MHD model of stars with 
complex magnetic fields (Long, Romanova & Lamb [20 IP] ) to 
the CTTS BP Tau. We investigate disc accretion onto the star 
by adopting the measured surface magnetic fields (D08) and 
other suggested properties of this star. The surface magnetic 
field is modeled as a superposition of a 1.2 kG dipole and 1.6 
kG octupole field, tilted by 20° and 10° with respect to the 
rotational axis and located at opposite phases (the phase dif- 
ference is 180°). We also take into account other parameters 
of BP Tau: its mass M* = O.7M ([Siess, Dufour & Foresting 
[2000] ), and radius R* = 1.95/2© ( [Gullbring et al.||1998p . Its 
age is about 1.5 Myr (D08), and its rotation period is 7.6 days 
( |Vrba et al.|19 86 ), which corresponds to a corotation radius 
of R cor ~ 7.5R*. The mass accretion rate derived from dif- 
ferent observations varies between M ~ 2.9 x 10~ 8 M Q yr~ 1 



also calculated the torque on the star and compared it with 
star's age. 

To understand the role of the octupole field in channeling 
the accreting matter, we compared our dipole plus octupole 
model of BP Tau with a similar model but with only the dipole 
component. 

Thus, we focus on: (1) 3D MHD modeling of accretion 
flows around CTTS BP Tau modeled with dipole plus octupole 
moments; (2) comparisons of accretion properties observed 
in simulations with observations of BP Tau, such as the shape 
and distribution of hot spots, mass accretion rates and more; 
(3) deviation of the simulated magnetic field from the poten- 
tial field. 

Section §2 briefly describes the numerical model used in 
this work. The simulation results are shown in §3. We end in 
§4 with our conclusions and some discussion. 



2 MODEL 



The global 3D MHD model originally developed by Koldoba et 
al. ([2002 ) and used for modeling stars with dipole fields ( |Ro 
man ova et al.||2003 2004a) was modified to include higher 
order components ( Long, Romanova & Lovelace 2007 , 2008 ; 



(e.g. [Gullbring et al.|[l998| and 9 x lQ-^Moyr -1 ( jSchmitt 



|et al.|2005| ). We performed a series of simulation runs at dif- 
ferent mass accretion rates in order to investigate the cases 
where the disk stops at different distances from the star. We 



Long, Romanova & Lamb 201 0[ |Romanova et al. 2010) to 
simulate disc accretion onto stars with complex fields. The 
MHD equations are solved in a reference frame co-rotating 
with the star. A viscous term is incorporated into the MHD 
equations (only in the disc) to control the rate of matter flow 
through the disk. We use the a— prescription for viscosity with 
a = 0.01. 

(i) Initial conditions. The simulation domain consists of a 
cold, dense disc and a hot, low-density corona, which are 
initially in rotational hydrodynamical equilibrium. The initial 
angular velocity in the disc is close to Keplerian. The angular 
velocity in the corona at any given cylindrical radius is set to 
be equal to that of the disk at that radius. 

(ii) Boundary conditions. At the inner boundary (the sur- 
face of the star), most of the variables A are set to have free 
boundary conditions, dA/dr — 0. The initial magnetic field 
on the surface of the star is taken to be a superposition of 
misaligned dipole and octupole fields. As the simulation pro- 
ceeds, we assume that the normal component of the field re- 
mains unchanged, i.e., the magnetic field is frozen into the 
surface of the star. At the outer boundary, free conditions are 
taken for all variables. In addition, matter is not permitted to 
flow into the region from the outer boundary. 

(iii) Simulation region and grid. We use the "cubed sphere" 
grid introduced by|Koldoba et al.| ( [2002 ) (see also Fig. 1 in 
Long, Roman ova & Lamb||2010| ). The resolution of the grid 
is 6 x 120 x 51 2 to simulate the accretion onto BP Tau in a 
simulation domain of 41i?*. We use a very high resolution 
near the star in order to resolve the complex structure of the 
octupolar component of the field. 

(iv) Magnetic field configuration. In our code we can model 
the magnetic field of the star by a superposition of three mul- 
tipole moments in (z = 1,2,3 for dipole, quadrupole and 
octupole respectively) which are tilted relative to the z— axis 
(which is aligned with the rotational axis ft) at different an- 
gles Si, and have different angles fa between the xz plane 
and ft — fii planes. For simplicity, fa is set to be 0, which 
means that the dipole moment \i\ is in xz plane. The general 
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Table 1. The reference values for CTTS BP Tau. The dimensional val- 
ues can be obtained by multiplying the dimensionless values from 
simulations by these reference values. Bo and the subsequent values 
below depend on pi . 



Reference Units 


pi = 1 


Ml =2 


Mi = 3 


M*(M ) 


0.7 








1.95 






Bi* (G) 


1200 






R (cm) 


3.9 x 10 11 






vo (cm s _1 ) 


1.5 x 10 7 






Po (days) 


1.83 






Bq (G) 


25.7 


12.9 


8.6 


po (g cm" 3 ) 


2.8 x 10- 12 


6.9 x 10- 13 


2.8 x 10- 13 


Mo (Moyr- 1 ) 


1.0 x 10~ 7 


2.6 x 10~ 8 


1.1 x 10~ 8 


Fo (erg cm _2 s _1 ) 


1.0 x 10 10 


2.6 x 10 9 


1.1 x 10 9 


N (g cm 2 s~ 2 ) 


3.9 x 10 37 


9.8 x 10 36 


4.6 x 10 36 



magnetic field configuration is discussed in greater detail in 
Long, Roman ova & Lamb] ( |2010| ) . The details of the magnetic 
field configuration in our BP Tau model are discussed in §3.1. 

(v) Reference units. The simulations are performed in di- 
mensionless variables A = A/A where A are reference val- 
ues. We choose the stellar mass M*, radius R* and the surface 
dipole field strength Bi* to build a set of reference values. 
The reference values are: length scale: R = i?*/0.35; ve- 
locity: vo = (GM*/R ) 1/2 ; time-scale: P = 2ttRo/v . The 
reference magnetic moments for dipole and octupole compo- 
nents are mi,o = B Rl and p 3 ,o = B Rl respectively, where 
Bo is the reference magnetic field. Hence, the dimensionless 
magnetic moments are: pi = pi/pi,o, P3 = p3/p3,o, where 
the dipole and octupole moments pi and p 3 of the star are 
fixed. We take one of the above relationships; for example, 
the one for the dipole; to obtain 



Bo 



pi,o 



05£i* 
Mi 



25.7 



£i* 
1.2kG 



G. (1) 



Hence, at fixed Bi*, the reference magnetic field depends on 
the dimensionless parameter pi . The reference density po and 
the mass accretion rate Mo also depend on this parameter: 

, 2 / i \ 2 

po — Bq/vo — 2.8 x 10" 



B lit 



Mo poVoRl 1.0 x 10" 7 



1.2kG 



B lir 
1.2kG 



cm° 

'Mo 
yr ' 



(2) 



(3) 



The dimensional accretion rate then is M = MM , where 
M is the dimensionless accretion rate.JDne can see that the 
dimensional accretion rate depends on M which we find from 
simulations, and the dimensionless parameter pi, which we 
use to vary the accretion rate. We change the dimension- 
less octupolar moment p3 in same proportion so as to keep 
the ratio pz/pn fixed. To find the ratio between the dimen- 
sionless moments, we use approximate formulae for aligned 
moments: pi = O.bBi+Rl and p 3 = 0.25B 3 *Rl (see Long, 
Romanova & Lamb 2010 ) and obtain for BP Tau the ratio 
Ji 3 /fa = B 3ie Rl/2B lie « 0.08 (where R* = 0.35). 

Other reference values are: angular momentum flux (a 
torque) No = poVqRo', energy flux E = p voRo', temper- 
ature T = IZpo/pQ, where 1Z is the gas constant; and the 



effective blackbody temperature T e ff,o = (po^o/ cr ) 1/4 > where 
a is the Stefan-Boltzmann constant. Tab. [T] shows the refer- 
ence values for CTTS V2129 Oph. In the subsequent sections, 
we show dimensionless values A for most of the variables and 
drop the tildes (~). However, we keep them in pi, JI3, and M 
because these are important parameters of the model. 

(vi) The magnetospheric radius. The truncation radius, r t , 
where the disc is truncated by the magnetosphere, could be 
estimated as (e.g., |Elsner & Lamb|1977| ): 

-1/7 



r t k(GM* 



'm- 2/7 /4 /7 , 



(4) 



where M is the accretion rate and mi is the dipole magnetic 
moment; k is a coefficient of order unity. For example, Long, 
|Romanova & L ovelace (2005) obtained k ~ 0.5 in numer- 
ical modeling of disc-accreting stars. For stars with known 
M*, R*, and dipole magnetic moment fi\, such as BP Tau, 
the accretion rate determines where the disc stops and how 
the matter flows onto the star. 



3 MODELING OF ACCRETION ONTO BP TAU 

We performed a number of simulation runs at different accre- 
tion rates and observed that the disk is truncated at different 
radii. We choose as our main case the one in which the disk 
stops sufficiently far away, but at the same time the accretion 
rate is not very low. Below, we describe this case in detail. We 
also show the results at higher accretion rates for comparisons 
(see §3.7). 



3.1 The modeled magnetic field 

D08 decomposed the observed surface magnetic field of BP 
Tau into spherical harmonics and found that the field is 
mainly poloidal with only 10% of the total magnetic energy 
in the toroidal field. The poloidal component can be approxi- 
mated by dipole (/ = 1) and octupole (/ = 3) moments with 
50% and 30% of the magnetic energy respectively. Other mul- 
tipoles (up to / < 10) have only 10% of the total magnetic 
energy. D08 concluded that the magnetic field of BP Tau is 
dominated by a 1.2 kG dipole and 1.6 kG octupole tilted by 
20° and 10°. The meridional angle between the ft — fj,± and 
ft — fi 3 planes is approximately 180°. In our model, we only 
consider the poloidal component. 

We convert the above parameters into dimensionless val- 
ues using our reference units and solve 3D MHD equations 
(see, e.g., |Koldoba et al.||2002| ), in dimensionless form. One 
of the important parameters of the model is pi which is used 
to vary the truncation radius in the dimensionless model (and 
the accretion rate in the dimensional model, see eq. 1). We 
find the second dimensionless parameter p 3 from the rela- 
tionship: p 3 /pi = £ 3 *#*/2£i* « 0.08 (where R* = 0.35). 
This ratio is fixed for fixed values of the dipole and octupole 
components, Bi* = 1.2kG and B 3 * = 1.6kG. From a num- 
ber of simulation runs at different pi, we choose pi = 3, 
p 3 = 0.24 to obtain a large enough magnetosphere as sug- 
gested by D08 and investigate this case in detail. Other pa- 
rameters of the magnetic field configuration of BP Tau are 
8i = 20°, 6 3 = 10° and 3 = 180°. 

Fig. [l] shows the components of the simulated magnetic 
field (B r , B0 and Be) in the polar projection down to the 
latitude —30°. One can see that the distribution of the radial 
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B(G) 



Figure 1. : Polar projections of the magnetic field components at 
the stellar surface in the dipole plus octupole model of BP Tau in 
spherical coordinates: radial magnetic field, B r ; azimuthal magnetic 
field, B^; and meridional magnetic field, Bq. The outer boundary, 
the bold circle and the two inner dashed circles represent the latitude 
of —30°, the equator, and the latitudes of 30° and 60° respectively. 
The red and blue regions represent positive and negative polarities of 
the magnetic field. 












- 






B(G)| 



Figure 2. The surface magnetic field in the dipole plus octupole 
model of BP Tau (fix = 3, fi 3 = 0.24, Bi = 20°, B 2 = 10°, 
03 = 180°) as seen from the equatorial plane (left panel), the north 
pole (middle panel) and the south pole (right panel) . The colors rep- 
resent different polarities and strengths of the magnetic field. 



component (left panel) is very similar to that obtained by D08 
in two observational epochs of Dec06 and Feb06 (see Fig. 14, 
left panels of D08) . In both cases, there is a strong positive 
pole at colatitudes of — 30° , a part of the negative octupolar 
belt at colatitudes of 60° — 90° and a part of the positive oc- 
tupolar belt at colatitudes of 90° - 120°. The distribution of 
the meridional component, Be, (right panel of Fig. [1} is qual- 
itatively similar to that of D08, though in the D08 plot the 
inner positive ring (red color in the plot) is weaker, while the 
negative ring (blue) is stronger compared with our model. 
The azimuthal component of the field is weak and shows a 
butterfly pattern of the negative and positive polarities in both 
the modeled and reconstructed fields. 

Fig. [5] shows three-dimensional views of the magnetic 
field distribution at the surface of the star. It can be seen 
that there are two antipodal polar regions of opposite polar- 
ity where the field is strongest. They approximately coincide 
with the octupolar high-latitude poles and their centers are 
located close to the octupolar axis fi 3 . Next to these regions, 
there are negative (blue) and positive (red) octupolar belts. 
Their shapes are more complex than those of the belts in pure 
octupole cases, where the belts are parallel to the magnetic 
equator (see |Long, Romanova & L amb 2010 ). This is because 
the dipole component strongly distorts the "background" oc- 
tupolar field. 

Although the octupole field is strongest at the surface of 
the star, it decreases more rapidly than the dipole field with 




Figure 4. 3D view of matter flow and the magnetic field distribution 
in the main case at t = 10. The left panel shows the distribution of 
the magnetic field lines. The middle panel shows one of the density 
levels, p = 0.84 x 10 -13 . The right panel shows the density distri- 
bution in the disc plane. The colors along the field lines represent 
different polarities and strength of the magnetic field. The thick cyan, 
white and orange lines represent the rotational axis and the dipole 
and octupole moments respectively. 



distance from the star. To investigate the role of the dipole 
and octupole components in channeling the accretion flow, 
we find the radius at which the dipole and octupole fields 
are equal. For this, we assume that both magnetic moments 
are aligned with the rotational axis and take the magnetic 
field in the magnetic equa torial planes: B\ — ni/r 3 an d 
B 3 = 3/i 3 /2r 5 (see Eqn. 1 in |Long, Romanova & Lamb|2010 ). 
Noting that the strengths of the field at the magnetic poles are 
B lic = and B 3 * = All 3 /RI, and equating B\ to B 3 , 

we find this radius: 



Teql 



3 £3* 

4 By* 



1/2 



R* 



(5) 



Substituting Bi* = 1.2kG and B 3 * — 1.6kG, we obtain r eq i ~ 
R*. Note that both dipole and octupole moments are tilted 
about the rotational axis, and hence the above formula gives 
only an approximate value for this radius. We suggest that 
this radius should be located slightly above the surface of the 
star to explain the dominance of the octupolar field seen in 
Fig.[l]and Fig.|2| 

Fig. [3] shows the initial magnetic field distribution near 
the star in our model. One can see the octupolar field com- 
ponent in the vicinity of the star. It also modifies the dipole 
field up to distances of r « 0.5R* (above the surface of the 
star) . The dipole field dominates in the rest of the simulation 
region. We should note that the dipole and octupole fields are 
equal at r eq i ~ R*, but the octupole component disturbs the 
dipole field up to larger distances. 



3.2 Matter flow and the magnetic field structure 

Here, we show results of 3D MHD simulations of matter flow 
onto the dipole plus octupole model of BP Tau with param- 
eters corresponding to the main case (fii — 3, fi 3 — 0.24) 
at time t — 10 when the system is in a quasi-stationary 
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state. Fig. [4] (middle panel) shows that the disk is truncated 
by the dipole component of the field, and matter flows to- 
wards the star in two ordered funnel streams. The right panel 
shows the density distribution in the equatorial plane. The 
left panel shows that the magnetosphere is disturbed by the 
disk-magnetosphere interaction. 

Fig. [5] shows slices of the density distribution and pro- 
jected field lines. Panel (b) shows the slice in the ^ — ft (or 
xz) plane. One can see more clearly that the dipole compo- 
nent of the field is responsible for disk truncation and matter 
is channeled in the ^ — ft plane. Panel (c) shows that the 
disk is stopped by the magnetosphere in the yz— plane. It is 
stopped at r w (6 — 7)R*. Panel (d) shows that matter flow is 
complex in the equatorial plane, and that matter comes closer 
to the star in the yz— plane. 

The gray line in Fig.[5]shows the magnetospheric radius, 
r m , the distance where the matter and magnetic stresses are 
equal: f3 = (p-\-pv 2 )/(B 2 /Sir) = 1. The magnetic stress domi- 
nates at r < r m ( [Ghosh & Lamb|1978[[l979||Long, Romanova| 
|& Lamb|20To] ). Fig. [5] shows that r m ~ which is smaller 
than the actual truncation radius r t ~ (6 — 7)R*. We should 
note that there are different criteria for the truncation radius 
and the above criterion usually gives the smallest truncation 
radius (see Bessola zet al.|20 08 for details). 

The external magnetic field evolves and deviates from 
the initial configuration due to disk-magnetosphere interac- 
tion. The left panels of Fig. [6] show the initial field distribu- 
tion at t = which corresponds to the potential field pro- 
duced by the star's internal currents. As the magnetic field 
evolves, the currents produced by the motion of plasma out- 
side the star become important and they change the poten- 
tial magnetic field significantly. The middle and right panels 
show the evolved field from different directions at t — 10. It 
can be seen that on a large scale the magnetic field lines wrap 
around the rotational axis of the star and form a magnetic 
tower (e.g. Lynden-Belf 1996[ |Kato, Hayashi & Matsumoto] 
|2004HRomanova et al.|2004b| ). This is a natural result of the 




p (9/cm d ) 



1.4E-14 2.6E-14 5.0E-14 9.4E-14 1.8E-13 3.4E-13 6.3E-13 1.2E-12 



Figure 5. Density distribution and selected field lines in different 
slices. Panels (a) and (b) show xz slices at t = and t = 10. Panels 
(c) and (d) show yz and xy slices at t = 10. The gray line shows the 
distance at which the matter stress equals the magnetic stress. 



magnetic coupling between the disk and star. The foot-points 
of the field lines on the star rotate faster than the foot-points 
threading the disk, which leads to differential rotation along 
the lines, and their stretching and inflation into the corona. 
The rotational energy is converted into the magnetic energy 
associated with these field lines. The evolved field structure 
on the large scale significantly differs from the potential field 
shown in the left panels. 

In the vicinity of the star the situation is different. The 
middle top panel of Fig. [6] shows that the magnetic field dis- 
tribution at t — 10 is almost identical to the initial field distri- 
bution. We conclude that the potential approximation is valid 
only within the parts of the magnetosphere where the mag- 
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Figure 6. Comparison of the initial (potential) magnetic field distribution at t = (left panels) with the field distribution at t = 10 (middle and 
right panels) . The top panels show the magnetic field in the vicinity of the star, while the bottom panels show the field distribution in the whole 
simulation region. The left and middle columns show the side view, while the right column shows the axial view of the field The color on the 
star's surface shows different polarities and strengths of the magnetic field. 




p (9/cm 3 ) 



1.3E-13 2.6E-13 5.3E-13 1.1E-12 2.2E-12 



Figure 7. The simulated accretion spots viewed from different di- 
rections at t = 10: from the equatorial plane (left-hand panel), the 
north pole (middle panel), and the south pole (right-hand panel). 
The color contours show the density distribution of the matter. The 
solid lines represent the magnetic moments of the dipole (/xi) and 
octupole (/X3). 
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Figure 8. The color background shows the distribution of the mag- 
netic field on the surface of the star. The blue contours show the en- 
ergy distribution in the accretion spot. The solid black lines show the 
dipole, octupole and rotational axes. 



netic stress dominates, that is, at r < 5R* in our case. The 
field distribution strongly departs from potential at larger dis- 
tances from the star. 

3.3 Accretion spots in BP Tau and in the Model 

D08 analyzed accretion spots from the observed brightness 
enhancements in chromospheric lines, such as Can IRT and 
Hei which presumably form in (or near) the shock front close 
to the stellar surface. The random flaring component, as well 
as the time-variable veiling component which reflects vari- 
ation of the intrinsic accretion rate were removed from the 
modeling. It was also suggested that matter flows in the vicin- 
ity of the strong magnetic field and hence Zeeman-splitting 
features were used for analysis. 

The corresponding maps of the local surface brightness 
are shown in Fig. 9 of D08 for two epochs of observations 



(Feb06 and Dec06). It can be seen that in both plots there is 
one bright spot located at colatitudes of 30° — 70° and cen- 
tered at C « 45°, while a lower-brightness area spreads up 
to a colatitude of 90°. In both cases the spots are elongated 
in the meridional direction. The Feb06 spot also has an an- 
tipodal spot of weaker brightness and of similar shape. The 
accretion filling-factor is shown in Fig. 13 of D08. It shows 
a spot located at colatitudes of 0° — 50° and centered at 
C « 10° - 20°. 

In our simulations the spots represent a slice of den- 
sity/energy taken across the funnel stream at the surface of 
the star. Hence, the spots show the distribution of these val- 
ues across the stream (Romanov a et al.|2004a| ). The physics 
of disk-magnetosphere interaction is expected to be more 
complex (e.g., |Koldoba et al.||2008| |Cranmer| [20081 12009| 
Bric khouse et al.|2010| ). However, this more complex physics 
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Figure 9. The simulated accretion spots in a polar projection down to 
colatitude 120° from the north pole. Left panel: density distribution; 
right panel: energy flux distribution. The equator is shown as a bold 
line. The dashed lines represent the latitude 30° and 60° respectively. 



strongly depends on the properties of funnel streams. Hence 
we call these spots "accretion spots" and show the density dis- 
tribution in spots and also the energy flux: F — pv • f[(v 2 — 
vl)/2 + 7^/(7 - l)p] distribution (right panel) at t = 10, 
where v* is the velocity of the star. 

Fig. [7] shows the density distribution at the surface of the 
star at t — 10. It can be seen that there are two antipodal 
spots which are centered in the \i\ — ft plane slightly be- 
low the dipole magnetic pole. The spots are elongated in the 
meridional direction (unlike the spots in a pure dipole config- 
uration, see also §3.4). 

Fig. [8] shows the density distribution in the spot (contour 
lines) overlaid on top of the magnetic field distribution. It can 
be seen that the spot is located near and below the dipole 
magnetic pole and hence its position is strongly influenced by 
the dipole component. The spot is located far away from the 
main, high-latitude magnetic pole, dominated by the octupo- 
lar component. We conclude that this low-latitude position of 
the spot results from the fact that the dipole field governs the 
matter flow, and only near the star, the octupolar component 
"comes to play" and the octupolar field redirects the spot's po- 
sition towards slightly higher latitudes, and changes the spot's 
shape. 

Fig. [9] shows a polar projection of the density distribu- 
tion (left panel) and energy flux. It can be seen that the ac- 
cretion spot is centered near the latitude of 40° and spreads 
between latitudes of 10° and 60°. These spots have a close 
resemblance to the spots observed by D08. In both cases, the 
spots are stretched in the meridional direction and are located 
between latitudes of 10° and 60°. Analysis of our spots at dif- 
ferent moments of time has shown that the positions of the 
spots vary only slightly with time (within 0.1 in phase). How- 
ever, the spots observed by D08 (their Fig. 9) are located at 
different phases. Feb06 observations show the main spot lo- 
cated at a higher phase compared with our spot, with a phase 
difference of 0.15-0.2. This is quite good agreement, consid- 
ering the fact the D08 method of the dipole moment phase 
reconstruction has an error in phase of 0.1-0.2 or larger. The 
phase difference is higher in the Dec06 observations where 
the simulated and observed spots are almost in anti-phase. 

To investigate the role of the dipole and octupole phases 
in the spot's position, we considered a few "extreme" cases, 
where the phase difference between the dipole and octupole 
moments is 4> — 0°, 4> — 90° and 4> — —90° (we have a 
phase difference of <j> — 180° in the main case). We observed 
that the accretion spot is similar in shape in all these cases, 
but is located at the phase corresponding to the phase of the 
dipole (in the \i\ — ft plane). We conclude that the phase of 
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Figure 11. Fraction of the star's surface covered by spots of density p 
and higher at t = 10. The solid and dashed lines represent the dipole 
plus octupole model and dipole model respectively. 



the accretion spots on the surface of the star gives a strong 
indication of the dipole component's position. 



3.4 Comparison with a pure dipole model 

To understand the role of the octupolar field component in 
our dipole plus octupole model of BP Tau, we investigate a 
similar model but with zero octupolar field, that is Bi* — 1.2 
kG, Bi = 20° and £ 3 * = 0. 

Fig.[To]shows xz slices of the accretion flow in the dipole 
plus octupole model (left panel) and pure dipole model (right 
panel) . One can see that the flow is wider in the dipole plus 
octupole model, because near the star the octupole field influ- 
ences the matter flow and redirects it in such a way that some 
matter flows towards the direction of the octupolar magnetic 
pole located at anti-phase with the dipole pole. The streams 
are narrower in the case of the pure dipole field. 

The bottom panels show that the accretion spots in the 
dipole plus octupole model are meridionally elongated, while 
in the dipole model they have their typical crescent shape and 
are elongated in the azimuthal direction (as usually seen in 
the pure dipole cases, e.g. Roman ova et al.| 2003 , 2004a]). 

We conclude that in BP Tau, the relatively weak octupolar 
magnetic field which dominates only very close to the surface 
of the star, strongly influences the shape and position of the 
accretion spots. This influence is not as dramatic as in another 
T Tauri star, V2129 Oph, which has a much weaker dipole 
component, and where the octupolar field splits the funnel 
stream into polar and octupolar belt flows (Roma nova et al.| 
|20T0l ). 



3.5 Area covered with spots 

D08 calculated the fraction of the star covered with visible 
chromospheric spots on BP Tau and concluded that they are 
spread over up to about 8% of the stellar surface, and cover 
about 2% of the surface, assuming one-fourth of each surface 
pixel is subject to accretion (see Figs. 9 & 13 in D08). In this 
paper we are interested only in the area of the spots spreading 
that is in the area covered by the funnel stream, and hence 
we take the larger value (8%). Figs. 9 and 13 of D08 show 
that the actual size of the spots depends on the brightness. 
Sometimes we take the brightest parts of spots (4-5)% and 
neglect the dimmer parts for comparison. In observations, the 
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Figure 10. Comparison of matter flow and accretion spots in the dipole plus octuple model of BP Tau (left panels) with a pure dipole case (right 
panels) at t = 10. The top panels show the density distribution (color background) and the sample poloidal field lines in the xz— plane. The 
bottom panels show the density distribution in the accretion spots as seen along the y— axis (panels a,c) and along the rotational axis (panels 
b,d). 




Figure 12. The density distribution of spots with a cutoff at different density levels p, and the corresponding spot coverage / for the sum of the 
visible and antipodal spots 



spot's coverage is calculated for one (visible) spot versus the 
whole area of the star, while in the simulations we take into 
account both spots, and hence the area in simulations is twice 
as large compared with observations. 

We observed from the simulations that the size of the 
accretion spot depends on the density (or energy) cutoff. We 
calculate the fraction of the star covered with the spots as 
f(p) = A(p)/4ttRI, where A(p) is the area covered with (all) 
spots of density p or higher. Fig. [TT] shows the distribution of 
f(p) versus p for the dipole plus octupole model and dipole 
model of BP Tau. It can be seen that for almost all chosen 
density levels, the accretion spots occupy a larger surface area 
in the dipole plus octupole model than in the dipole model. 
This result is in good agreement with that obtained from Fig. 
M 

We choose several cutoff densities and show the spots 
in Fig. [12] One can see that in both simulations and obser- 
vations, spots cover different areas depending on the bright- 
ness/energy flux levels ranging from 14% up to 2% in simula- 
tions (for two spots) . In both observations (see Fig. 9 of D08) 



and simulations (see our Fig.[To]), the brightest parts of spots 
are located at 30° < 6 C < 60°. 

It was predicted that the area covered by spots in the case 
of complex magnetic fields should be smaller than in the pure 
dipole c ase (e.g.,|Mohanty & Shu|2008[ |Gregory et aL|2 008, 
see also Calvet & Gullbring 1998 }. This paper and our previ- 
ous work (Long, Romanova & Lovelace 2008 ) show that the 
area / could show a complex trend for mixed dipole and mul- 
tipole configurations, depending on how the multipole field 
redirects the accretion flow. 



3.6 Accretion rate in BP Tau and in our model 

The mass accretion rate of BP Tau obtained from observations 
is not uniquely determined, and different authors give differ- 
ent results. For example, Gullbring et al. (1998 ) estimated the 
mass accretion rate to be 2.9 x lO _8 M yr _1 using the lumi- 
nosity of L ~ {GM*M/R*)(l - R*/R in ), where R in is the 
inner radius of the disk. Other estimates of the mass accre- 



tion rate include 1.6 x 10~ 7 M©yr _1 (Valenti & Johns-Krull 
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2004), 9 x lO-^Moyr" 1 flSchmitt et al.|2005| ). |Calvet et al. 
( 2004| ) compared the mass accretion rates of intermediate- 
mass, (1.5 - 4)M , and low-mass, (0.1 - l.O)M , T Tauri 
stars and concluded that the average mass accretion rate for 
intermediate-mass T Tauri stars is about 3.0 x lO _8 M yr _1 , 
while for low-mass CTTSs (like BP Tau) it is about 5 times 
smaller, that is, 6.0 x lO~ 9 M yr _1 . 

We obtained from simulations the dimensionless accre- 
tion rate of M « 0.13. Using the reference value for M from 
Tab.[l] M = 1-1 x lO^Moyr -1 (calculated for Jli =^3), 

we obtain the dimensional accretion rate: M — M M w 
1.4 x lO _9 M yr _1 . This value is higher than the accretion 
rate suggested by Sch mitt et al.| ( |2005| ) but lower than the 
other estimates. Given that BP Tau has the low mass of 
- O.7M (D08) (or O.49M [Gullbring et al.|1998p , the mass 
accretion rate from our simulation models may still be in the 
reasonable range according to the |Calvet et al.| ( [2004| ) analy- 
sis and discussion. 

Here, we should note that the "main model" considered 
above is one of the suggested models, where the disk is trun- 
cated at large distances. Below, we discuss additional models, 
where the disk is truncated at smaller distances and the ac- 
cretion rate is higher. 



3.7 Modeling of BP Tau at higher accretion rates 

We performed two additional simulation runs at higher ac- 
cretion rates (at smaller values of parameter Jli) keeping 
the ratio between the dipole and octupole moments fixed, 
Jl 3 /Jli = 0.08. The parameters of these new models are the 
following: Jn = 2, JI3 = 0.16 and Jli = 1, JI3 = 0.08. 

Fig. [13] (middle and right panels) shows the matter flow 
and accretion spots for these two cases, while the left panel 
shows results for the main model (Jli = 3,^3 = 0.24) for 
comparison. One can see that in the new models, the disk is 
truncated at smaller distances from the star: r t ~ (5 — 6)R* 
(for Jli = 2) and r t ~ 3.6R* (for Jli = 1). The bottom panels 
of Fig. [13] show that when the disk comes closer to the star, 
the spots are still at high latitudes, though they become longer 
in the meridional direction. This is probably because the oc- 
tupolar component has a stronger influence on the spots: the 
spots move closer to the octupolar magnetic pole. Note that 
no octupolar ring spots appear (like in V2129 Oph, see R10). 
This is probably because the disk truncation radius is still far 
away from the area where octupole has a strong influence on 
the flow, which is closer to the star. It is interesting to note 
that in the pure dipole case, the spots move towards lower 
latitudes when the disk comes closer to the star. Here, we see 
the opposite: parts of the spots move to higher latitudes, be- 
cause the octupolar component becomes more significant at 
smaller truncation radii. 

The dimensionless massjaccretion rate which we ob- 
tain from our simulations is M w 0.075 (for Jli — 2) and 

M « 0.085 (for Jli = 1). We take reference values Mo for dif- 
ferent Jli from Tab. [T] and obtain corresponding dimensional 
accretion rates M (see Tab. [2]) which are 2.5 x lO~ 9 M yr _1 
and 8.5 x lO _9 M yr _1 respectively. We see that the accretion 
rate for Jli = 1 is high and close to many of the observed 
values. However, the disk comes too close to the star and the 
truncation radius is much smaller than the corotation radius 
(Rcor ~ 7.5i?*), which would mean that the star is not in 
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Figure 13. Simulations at different accretion rates at t = 8. (a) 
main case, Jli = 3,/i3 = 0.24; (b) Ji\ = 2,^3 = 0.16; (c), 
Jli = 1,Jls = 0.08. The top panels show the density distribution 
(color background), the poloidal field lines (red lines) and the (3 = 1 
line (gray) in the xz— plane. The bottom panels show the energy flux 
distribution in the accretion spots in a polar projection. 



the rotational equilibrium state, which is unlikely The model 
with Jli = 2 is similar to the main model: the disk is truncated 
far away and the accretion rate is twice as high, but still it is 
lower than many values given by observations. This model is 
slightly better than the "main" model. 

In the case of a pure dipole, we obtain from simulations 

M « 0.12 and M ^ 1.3 x 10 " 9 M yr _1 . [2] summarizes the 
results obtained in our simulation models and some observa- 
tional properties of BP Tau for comparison. 



3.8 Angular momentum transport 

We calculated torque on the surface of the star for the above 
three models. The torque associated with the magnetic field 
Nf is about 10-30 times larger than the torque associated 
with matter and thus dominates (see also |Romanova et al.| 
|2002||2004al|Long, Romanova & Lovelace|20071 12008| ) . The 
dimensionless torque obtained from simulations, N , varies in 
time around some average value. Tab. [2] shows these values 
and values corresponding to these variations. We observed in 
simulations that the torque is negative in most of the models 
because the disk is truncated at the distances comparable with 
the corotation radius and therefore the star spins-down. In a 
model with Jli = 1, the disk is truncated at the distance of ~ 
0.5R C or and it spins the star up. To estimate the dimensional 
torque, Nf, we take the largest absolute value of Nf from 
Tab. [2] and take into account the reference values of iVo from 
Tab.[l] Tab.|5]shows the dimensional torque. 

We estimate the time-scale of spinning-down of BP Tau. 
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Table 2. Comparison of observational and modeled results shown in D08 and results from our simulations. 



Models 




n 


M 


MCMoyr- 1 ) 


N f 


Nf (g cm 2 s~ 2 ) 


r (yr) 




D08 & other measurements 


>4R* 




1CT 9 - 1CT 7 










Model (pi = 3, JI3 


= 0.24) 


(6 - 7)#* 


0.13 


1.5 x 10- 9 


-0.03 ±0.014 


-1.9 x 10 35 


1.6 x 10 7 


(spin-down) 


Model (pi = 2, p 3 


= 0.16) 


(5.4 - 6.4)fl* 


0.075 


1.9 x 1CT 9 


-0.005 ±0.01 


-1.5 x 10 35 


2.1 x 10 7 


(spin-down) 


Model (pi = 1, JI3 


= 0.08) 


3.6#* 


0.075 


7.7 x 10- 9 


±0.002 ± 0.004 


+ 2.3 x 10 35 


1.1 x 10 7 


(spin-up) 


Model (pi = 3, ps 


= 0) 


(6 - 7)fl* 


0.12 


1.4 x 1CT 9 


-0.045 ±0.015 


-2.4 x 10 35 


1.3 x 10 7 


(spin-down) 



The star's angular velocity is (at period P = 7.6 days) 
Q = 2n/P « 9.6 x 10 _6 s _1 , its angular momentum is 
J = kM+RlQ. = 2.5 x 10 50 k gcm 2 /s, where k < 1 (we 
take k = 0.4 for estimations). The spin-down time-scale, 
r = J/Nf, was estimated for different models and is shown in 
Tab. [2] One can see that the spin-up/down time-scale is about 
an order magnitude larger than the age of BP Tau (1.5 x 10 6 
years, see references in D08). Hence, the torque obtained in 
simulations is not sufficient to regulate the spin of the star 
at the considered epoch. It is possible that a star lost most of 
its angular momentum at the earlier stages of its evolution. 
Note that similar low-torque situation has been observed in 
modeling of another CTTS, V2129 Oph. 



4 CONCLUSIONS AND DISCUSSIONS 

We performed global 3D MHD simulations of accretion onto 
a model star with parameters close to those of the classical T 
Tauri star BP Tau, and with the magnetic field approximated 
with a superposition of slightly tilted dipole and octupole mo- 
ments with polar magnetic fields of 1.2kG and 1.6kG which 
are in anti-phase (D08) . In this star the dipole field dominates 
almost up to the surface of the star and determines the ma- 
jority of the observational properties. 

We performed a number of simulation runs for different 
truncation radii of the disk, and chose one of them where the 
disk is truncated at r « (6 — 7)R* which is sufficiently close 
to the corotation radius R cor ~ 7.5.R* and investigated this 
case in greater detail. We observed that the dipole component 
of the field truncates the disk, and matter flows in two funnel 
streams towards the dipolar magnetic poles. However, near 
the star the flow is slightly redirected by the octupolar com- 
ponent towards higher latitudes, and this affects the shape 
and position of the accretion spots: the spots are stretched 
in the meridional direction and are centered at higher lati- 
tudes compared with spots in the pure dipole case, which are 
latitudinally-elongated (crescent-shaped) and are centered at 
lower latitudes. 

The spots are located near the hi — ft plane, where both 
the dipole and octupole moments are situated, and are in anti- 
phase. Experiments with different relative phases between the 
dipole and octupole moments have shown that the spots are 
always located in the meridional plane of the dipole moment. 

The spot's position slightly varies in phase with the ac- 
cretion rate, but the variation is small, about 8° at the most. 
Note that in the pure dipole case, the spot may rotate about 
the magnetic pole (and hence may strongly change the phase) 
if the dipole is only slightly tilted about the rotational axis, 
Bi < 5° ( Romanova et al.|2003|[200^|Bachetti et al.|2010) . 
However, in BP Tau, the tilt is larger, Bi « 20° and an oc- 



tupolar component is present and is significant near the star. 
Both factors lead to the restriction of such a motion. 

The meridional position and shape of spots observed 
in simulations is similar to those observed by |Donati et al.| 
( 2008 ) (see their Fig. 9 for brightness distribution in spots re- 
constructed for the Dec06 or Feb06 epochs). However, they 
are at different phases compared with our spots. This may be 
connected with the relatively low accuracy of the phase recon- 
struction of the dipole component from the surface magnetic 
field, or some other reason. 

The accretion rates obtained in different models are in 
the range of M w (1.5 - 7.7) x lO" 9 M yr _1 . It is lower than 
most of the M values derived from observations which vary 
in the range of 9 x 10~ 10 and 1.6 x 10~ 7 solar mass per year 
and depend on the approaches used for the derivation. The 
smallest accretion rate obtained in simulations corresponds 
to the case where the disk is truncated at r t ~ R cor , and a 
star slightly spins down, while the largest to the case where 
r t ~ 0.5 Rcor and a star spins up. We should note that at 
larger accretion rates, say at 1.6 x 10~ 7 solar mass per year, 
the disk will come very close to the star, and the torque would 
be much stronger. Such a state seems unlikely. Hence, if the 
accretion rate is very high, then we should suggest that the 
dipole component should be larger than that derived by D08. 

The torque obtained in simulations is small and the 
time-scale of the spinning-up/down is an order of magni- 
tude smaller compared with the age of BP Tau. This torque 
is not sufficient to support a star in the rotational equilib- 
rium state, where a star spins up or spins down depend- 
ing on the accretion rate but has a zero torque on average 
(e.g.,|Ghosh & Lamb|1979||Konigl|1991||Cameron & Campbell| 



1993 Long, Romanova & Lovel ace 2005). We assume that 
this small torque matches a small accretion rate obtained in 
simulations because the torque generally correlates with the 
accretion rate (e.g., Romano va et al.|2002| ). 

We suggest that a star may lose the majority of its an- 
gular momentum at earlier stages of its evolution due to the 
"propeller" effect (e.g., Romanova et al. 2005; Ustyugova et 
al. 2006), stellar winds ( [Matt & Pudritz||2005| 2008), or by 
some other mechanism. 

Earlier, we performed global 3D MHD simulations of ac- 
cretion onto V2129 Oph ( Roma nova et al.|2010| ) which have 
shown that in the case of a strong octupolar component, parts 
of the octupolar belt spots can be visible and can dominate at 
sufficiently high accretion rates. This is not the case in BP Tau, 
where the octupolar component dominates only in the close 
vicinity of the star. 

Disk-magnetosphere interaction leads to inflation of the 
external field lines and formation of a magnetic tower. Sim- 
ulations show that the potential approximation used in ex- 
trapolation of the magnetic field from the surface of the star 
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to larger distances (e.g. Jardine et al. 2006, D08, Gregory 
|et al.|20io| ) is valid only inside the magnetospheric (Alfven) 
surface, where the magnetic stress dominates. At larger dis- 
tances, the magnetic field distribution strongly departs from 
potential. 
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